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Abstract. Quantitative estimates of lensing probabilities must be self- 
consistent. In particular, for asymptotically isothermal models: (1) using 
the (3/2)^/^ correction for the velocity dispersion overestimates the ex- 
pected number of lenses by 150% and their average separations by 50%, 
thereby introducing large cosmological errors; (2) when a core radius is 
added to the SIS model, the velocity dispersion must be increased; and (3) 
cross sections and magnification bias cannot be separated when computing 
the lensing probability. When we self-consistently calculate the effects of 
finite core radii in fiat cosmological models, we find that the cosmological 
limits are independent of the core radius. 



1. Introduction 

Asymptotically isothermal potentials are consistent with most data on 
gravitational lenses. They explain the observed numbers of lenses (Maoz 
& Rix 1993, Kochanek 1993), fit most observed image configurations (e.g. 
Kochanek 1991a), and are consistent with stellar dynamics (Kochanek 1994, 
Franx 1993). Photometry of early type galaxies (Tremaine et al. 1994), and 
the absence of central images in most lenses (Wallington & Narayan 1993, 
Kassiola & Kovner 1993) suggest that the lens potentials have a small or 
vanishing core radius. There is, however, a persistent myth (generated in 
part by the author) that small core radii can dramatically alter the ex- 
pected number of lenses without other observational consequences. We can 
trace the current versions of this myth to inconsistencies in either the dy- 
namical normalization of the models or the calculation of the magnification 
bias. In this technical memorandum we briefiy explore these consistency re- 
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quirements and the resulting effects of a finite core radius on cosmological 
limits. 

We confine ourselves to a circular isothermal density distribution with 
p = ajjj^/2'iTG(r^ + s^) where (Jdm is the velocity dispersion of the dark 
matter, and s is the core radius (Hinshaw & Krauss f987). The lens deflects 
rays by 

dc^ _ ^ (r^ + s^y/^ - s 
dr r 

where b = ^i^{ctdm I c)^ D^g I Dqs-, Dls and Dqs are proper motion or an- 
gular diameter distances between the lens and the source and the observer 
and the source respectively, and (j) is the two-dimensional lensing potential. 
(The alternate softened isothermal model with (]) = b(r^ + s^)^/^ intro- 
duced by Blandford & Kochanek (f987) is almost indistinguishable from 
the Hinshaw & Krauss (1987) model if its core radius is twice as large.) 

The lens is supercritical and able to generate multiple images if /3 = 
s/b < 1/2, with a tangential critical line at r_|_ = 5(1 — 2/3)^/^ and a radial 
critical line at r_ = b[(3 - (3^/2 - (3^^^(4 + /3)^/2]i/2_ r^^^ caustics are at 

= and M_ where the cross section is r = ttm^ = 7r5^[l -|- 5/3 — /3^/2 — 
^1/2(^4 _l_ ^^3/2^2] (Hinshaw & Krauss 1987). The cross section declines very 
rapidly with /3 and near the threshold of /3 = 1/2 — e it declines as r oc e"^ 
if b is held fixed. If we assume a constant comoving core radius the cross 
section can be integrated analytically to compute the optical depth (Krauss 
& White 1992). 

2. Dynamical Normalizations 

The first question we must address is the normalization of the singular 
model (s 0). Historically, Turner, Ostriker, & Gott (1984) argued that 
if the central velocity dispersion of the stars is ac then the dark matter 
should have velocity dispersion (T_dm = (3/2)^/^(Tc. However, Franx (1993), 
Kochanek (1993, 1994), and Breimer and Sanders (1993) show convincingly 
that real galaxies do not satisfy the assumptions used by Turner, Ostriker, 
& Gott (1984), and that for real galaxies (Jdm — <^c- Kochanek (1994) fit a 
sample of 37 early type galaxies from van der Marel (1991) and found that 
the best fit estimate was (Jdm* = 225 ± 10 km s~^ for an galaxy. 

All existing studies of the effects of a core radius on lens statistics have 
added a core radius while leaving the velocity dispersion (Jdm or b un- 
changed. It is clear, however, that s and (Jdm must be correlated. Adding a 
core radius reduces the mass near the center of the galaxy, and the velocity 
dispersion must increase compared to its value in a singular model to main- 
tain either the stellar velocity dispersions or the average image separations 
fixed. As a model calculation, we compute the average line-of-sight veloc- 



A TECHNICAL MEMORANDUM ON CORE RADII 



3 



ity dispersion inside one effective radius Re assuming a Hernquist (f990) 
distribution (z/ oc r~^(r + «)~"^) for tfie stars witfi a ~ 0.45i?e- Witfi tfie 
assumption tfiat tfie vefocity dispersion tensor of tfie stars is isotropic, tfie 
dark matter dispersion increases as (Jdm oc f + 2(s/i?e) witfi tfie addition 
of a core radius. 

For a more reafistic modef we use van der Maref (f99f) sampfe and fit 
isotropic dynamicaf modefs to eacfi gafaxy, assuming tfie core radius is a 
constant fraction of tfie estimated effective radius for eacfi gafaxy, and tfiat 
tfie vefocity dispersion scafes as L/L^ = {ctdm I '^DM*^ ■ Tfie surface 
of tfie fit to tfie observed vefocity dispersion profifes is sfiown in Figure f , 
and tfie dasfied fine is tfie scafing faw estimated from tfie Hernquist (f990) 
modef. Modefs witfi farge core radii cannot fit tfie data because of tfie 
contradiction between a fiomogeneous core and a steepfy rising fuminos- 
ity profife. Tfie formaf 95% confidence upper fimit on tfie core radius is 
s* ^ 0.08i?e or ^ 300/i~^ pc for i?e* = 4/i~^ kpc. For a core radius of 
= f 00/i~^ pc tfie fractionaf increase in tfie vefocity dispersion is 7.5% or 
f 7 km s~^, fess tfian tfie uncertainty in tfie vafue of udm*- Nonetfiefess, its 
effects on modefs witfi a finite core radius are striking; it produces a 33% 
increase in tfie expected number of fenses if we keep tfie ratio of tfie core 
radius to tfie criticaf radius fixed constant). 

Seff-consistency in fens modefs afso requires a vefocity dispersion tfiat 
increases as tfie core radius becomes farger since tfie average image separa- 
tion must stay fixed as tfie core radius increases. Tfie image separation is 
approximatefy twice tfie tangentiaf criticaf radius of tfie fens (A^ ~ 2r_|_ = 
25(f — 2/3)), so tfiat if tfie core radius increases (farger /3), tfie onfy way to 
maintain constant average image separations is to afso increase tfie average 
vefocity dispersion (farger h). If we modef tfiis by keeping tfie tangentiaf 
criticaf radius r_|_ = 5(f — 2/3) fixed, tfien tfie fens cross section r decreases 
as r oc instead of . 

3. Magnification Bias 

Seff-consistent cafcufations of tfie fensing probabifity sucfi as Kocfianek & 
Bfandford (f987), Kocfianek (f99fb, f993), Waffington k Narayan (f993), 
Kassiofa & Kovner (f993), and (in most regimes) Maoz & Rix (f993) au- 
tomaticaffy incfude tfie effects of tfie core radius on tfie magnification bias, 
but most treatments of softened isotfiermaf modefs examined onfy tfie ef- 
fects of a core radius on fensing cross sections (e.g. Dyer f 984, Bfandford & 
Kocfianek f987, Hinsfiaw k Krauss f987, Krauss k Wfiite f992, Fukugita 
k Turner f 99f , Fukugita et af. f 992, Bfoomfiefd- Torres k Waga f 995). Core 
radii fiave a powerfuf effect on tfie cross section for muftipfe imaging. How- 
ever, using tfie cfiange in tfie cross section grossfy overestimates tfie effects 
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Figure 1. (Left) Contours of the for the dynamical fits. The light solid lines show the 
68%, 90%, 95%, and 99% confidence limit changes on Ax^ for the fit to the sample. The 
dashed line shows the expected scaling of <tbm* with s/ Re if we keep the average velocity 
dispersion interior to Re fixed in the Hernquist /softened isothermal sphere dynamical 
model. 

Figure 2. (Right) Variation in cross section (dashed line) and lensing probability (solid 
line) with core radius s for a lens with <jdm = 250 km . The values are normalized 
to unity at the minimum core radius. The heavy solid line shows the lensing probability 
excluding image systems with detectable central images using the same normalization 
as for the total probability. The results are given for the average over the quasar data 
sample including selection effects. 

of a finite core radius on the lensing probability in bright quasar samples. 
The core radius first eliminates images with low total magnifications, but 
the bright quasar samples are dominated by highly magnified images of 
fainter quasars and magnification bias significantly reduces the effects of 
adding a core radius on the probability. 

We can understand this analytically in the Hinshaw & Krauss (f987) 
model. Inconsistent models of the effects of a core radius estimate the lens- 
ing probability by using the optical depth multiplied by the magnification 
bias for the singular model. When two images are merging on the radial 
caustic, the third image is located at Tout = '2(3u-/r'^, and the average 
magnification produced by the lens is (M) = t^^^/m^ = 4/3^ /rl. When the 
core radius is small, the average magnification is 4, but near the thresh- 
old the average magnification diverges as (M) oc e~^. The magnification 
probability distribution is approximately P(> M) = ((M)/2M)^ when 
M > (M)/2 for fold caustic statistics (e.g. Blandford & Narayan f986). 
If we assume a single power law quasar number counts distribution with 
dN/dm oc 10"'(™-™«) then the magnification bias varies with the average 
magnification as B{m) oc (M)o^"' for 7 < 0.8. As the core shrinks, the av- 
erage magnification increases, which drives up the magnification bias. The 
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Figure 3. Cosmological effects of finite core radii. Tfie ieft panei sfiows contours as 
a function of tfie ratio of tlie core radius to tlie effective radius and tfie cosmoiogicai 
modei. Tfie iigfit soiid fines are tfie constraints from fensing afone, tfie dasfied fines 
adds tfie vefocity dispersion prior as a function of core radius estimated in §2, and tfie 
fieavy sofid fines add tfie goodness of fit to tfie dynamicaf data. Tfie rigfit panef sfiows 
dependence of ctdm* on tfie core radius in an fio = f cosmofogy. Tfie figfit contours use 
onfy tfie fens data, and tfie dasfied fines incfude tfie prior probabifity distribution for 
o'DM* estimated from dynamicaf modefs. Tfie contours are drawn at tfie 68%, 90%, 95%, 
and 99% confidence fevefs for one parameter. 

lensing probability, including the change in the magnification bias, varies 
as TB{m) (X e'^~^^ not r oc e"^. For large average magnifications the effective 
value of 7 is the faint slope of the quasar number counts, 7 ~ 0.27±0.07, and 
TB(m) oc fi-65±o.35_ bright quasars the increase in the bias is greater 
because of the steeper number counts slope, and the effects of the core 
radius are still smaller. 

To emphasize this point. Figure 2 shows the relative variation of the 
cross section and the true lensing probability including magnification bias 
for a lens with (Jdm = 250 km s~^ as a function of the core radius averaged 
over the full quasar data sample. For a core radius of s = 100/i~^ pc using 
the cross section instead of the true probability underestimates the lensing 
probability by about 40%. This comparison overestimates the effect of a 
finite core radius because we did not include the dependence of a dm on s. 

4. The Cosmological Effects of Softened Isothermal Spheres 

We assume that the core radii of galaxies are proportional to their effective 
radii s = and the models are characterized by a fixed ratio of 

s/Re- To simplify the calculations, we set the "Faber- Jackson" exponent 
to be 7 = 4, so the core radius varies with the velocity dispersion as s = 
s*icrDM / '^DM*Y'^ = s^{L I L^y-'"^ (see Kochanek 1995 for details). 
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The right panel of Figure 3 shows the best fit value of (Jdm* as a function 
of in an i7o = 1 cosmology. As expected from §2, the velocity disper- 
sion increases as the core radius increases with (Jdm* oc 1 + s/Rg. This is 
shallower than the slope seen in the dynamical models. The left panel of 
Figure 3 shows the dependence of the cosmological limits for a fiat universe 
(fio + Ao = 1) on the core radius using only the lens data, the lens data 
combined with the prior probability distribution for (Jdm* derived from the 
dynamical model of §2, and finally the lens data, the dynamical velocity 
dispersion prior, and the likelihood of the dynamical model. If we use only 
the lensing data, the cosmological limits are nearly independent of the core 
radius - this is a radically different picture of the effect of a core radius 
than that found in inconsistent calculations. 
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